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A universal method of strictly calculating self-consistent fields of realistic 
plasma particles could be strictly derived from three basic tools in theoretical 
plasma physics: particle simulation, Vlasov-Maxwell theory and fluid theory. 
PACS: 52.65.-y, 52.35.-g. 
Plasma physics is a physical branch about many charged particles interacting through 
their self-consistent fields. In its earlier developing stage (about 1940s~ 1960s), many the- 
oretical methods [1-6] which are successful in other elder physical branch such as neutral 
gas physics and fluid mechanics were transplanted into this younger branch and rapidly 
built up the basis of this new branch. However, almost no one doubts whether these trans- 
planted methods are appropriate for plasmas where numerous charged particles are corre- 
lated through their self-consistent fields. More important, in above-mentioned transplanted 
methods the plasma self-consistent fields is never strictly calculated but is indeed treated 
by various (obvious and hidden) approximations. 

Although people have realized the importance of strictly calculating plasma self- 
consistent fields, some realistic adverse factors prevents this goal being thoroughly achieved. 
Let us comment three basic tools in plasma physics one by one. 

a) particle simulation. [7-9] 

This basic tool is to solve 2A^ + 4 equations describing N realistic particles (or macropar- 
ticles) 
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d t E(R, t) = VxB(R J t) + J2 dtn (t) 6 (r 4 (t) - R) ; (P.l) 

i 

d t B(R, t) = -V x E(R, t); (P.2) 

V-E(R,t) = J2H^(t)-R); (P.3) 

i 

V-B{R,t) = 0; (PA) 



d t = = S(r< (t) , t) + dtrj (t) x 5(r, (t) , t); (P2i+4) 

1 - (t)] 2 

^ = (P2i+5) 



It is in principle a strict method. But in practice, because N is nearly an astronomical 
figure, people often resort to an approximation method, which is often called Particle-In- 
Cell (PIC) method. The PIC method is to approximate the solution of above 2N + 4 
equations with that of 2N / R merge + 4 equations, where R merge > 1 measure how many 
realistic particles are merged into a so-called macroparticle, and to alternatively updating 
N/R merge macroparticles' information (i.e., position and velocity) and (E,B). 

b) Vlasov-Maxwell theory. [1-6] 

This basic tool is to solve 5 equations 

d t E(R, t) = V x B(R,t) + J vfd 3 p; (V.l) 

d t B(R, t) = -Vx E(R, t); (V.2) 

V-E(R,t) = J fd 3 p-N z ; (V.3) 

V-B(R,t) = 0; (V.4) 

d t f + vVf-[E + vxB}-d p f = 0. (V.5) 

It is also in principle a strict method. But in practice, because / is defined over a 6-D phase 
space and hence corresponds to too huge data mount, if above 5 equations are solved by 
alternatively updating / and (E,B), updating / will be very time-consuming. This basic 



tool is therefore less applied than the approximation version of particle simulation, i.e., PIC 
method. 

c ) fluid theory. [1-6] 

This basic tool is to solve 5 + 1 equations 

d t E(R, t) = V x B (R,t) + nu fl - 
dtB(R,t) = -V x E(R,t); 

V ■ E(R, t) = n-Ni] 

V • B(R, t) =0; 

d t pfi +Ufi ■ Vp/; = E + Ufi x B + thermal pressure /density; 

Assumed ther mod ynamical state equation (about thermal pressure/density) 

Because of assumed thermodynamical state equation, this basic tool is often viewed as 
inferior than other two tools. The self-consistent fields (E, B) obtained from this tool is 
hence taken as less reliable than its counterparts obtained from other tools. 

Following text will show in details that if above three basic tools are in their respective 
strict forms, they will agree with each other to yield a strict method of calculating (E, B) 
of realistic particles. 

a) for particle simulation 

We could rewrite any relativistic Newton equation, for example Eq.(P.2i+4), as 

dkn (t) 



= d, 



d t r l it) 



E( ri (t),t)-dtri (t) x B(n (t),t) 



-d. 



u (r; it) , t) 



+ 



d 



1 - [dtTi {t)f y/l-Hr, (t),t)] 

u (n (t) , t) 



- [dtn {t)-u{n (t),t)] x B( ri (t),t) 



y/l-[u( ri (t),t)f 



E( n (t),t)-u(ri (t),t) x B( ri (t),t) 



(fp.l) 



3 



and note a fact that it is valid for arbitrary value of d t r j (t) , or arbitrary value of A = d t r , (t) — 
u (ji (t) ,t). Eq.(fp.l) is of a binary-function type general form: = function(varl,var2), 
where varl and var2 are independent variables. Thus, timing Dirac function 8 (A) at both 
side of Eq.(fp.l) and then integrating over A, we could obtain 







[right-side terms of Eq.(fp.l)] * 6 (A) dA 



d f 



u (r; (t) , t) 



- E (n (t) ,t)-u (n (t) , t) x B{n (t) , t) 



(fp.2.b) 



(t) } t)} z 

Here, we have utilized, when deducing Eq.(fp.2.b), following property of the Dirac function: 
x8 (x) = 0. This property immediately leads to [dtx] * 6 (x) = —x * [d t 6 (x)]. Noting the 
property of d t S (x): d t S (x) = if d t x = 0;and d t S (x) = 8 (x) if d t x ^ 0, (i.e., if x varies 
with respect to t, x-value will derivate from and corresponding 8 (x)-value will also jump 
from oo to 0), we could find that there are x * [d t 8 (x)] = if d t x = and x * [d t 8 (x)] = 
x8 (x) = if d t x ^ 0, i.e, no matter what dfX-value is, there is always x * [d t 8 (x)] = 0, 
and hence [dtx] * 8 (x) = 0. The integral J [right-side terms of Eq.(fp.l)] * 8 (A) dA includes 
terms of a general form J d t A * 8 (A) dA. These properties of the Dirac function lead to 
/ d t A * 8 (A) dA = 0. 

Subtracting Eq.(fp.l) and Eq.(fp.2.b), we have 







dt. 



d t Ti (t) 



d 



U {Ti (t) , t) 



V 1 - Wi it)f y/l-[u( ri (t),t)f 



[dtn {t)-u{n (t),t)] x B{n (t),t). 

(fp.2.a) 



Therefore, any solution of Eqs.(P) is also that of following equation set of 2N + 5 members 



u f i {n (t) , t) 



[Ufl {Ti (t),t)f 



- E (n (t) , t) - u f i {n (t) , t) x B{n (t) , t) 



(2.0) 



d t E(R, t) = VxB(R,t) + J2 d t r * (*) 5 in ^) ~ R ) ; (2- 1 ) 

i 

d t B(R, t) = -V x E(R, t); (2.2) 



V-E(R,t) = Y / 6{r i (t)-R); (2-3) 

i 

V-B(R,t) = 0; (2.4) 



, d t n(t) Ufi(ri(t),t) 
at — , — Ctf 



^l-[dtn(t)} 2 y/i-[u fl (n (t),t)] 2 _ 



[dtU (t) -Ufi {n (t) , t)} x B{n (t) , t) (2.2i+4) 
Vi = dtn (2.2i+5) 



where (R, t) = J2 i€r .( t )=R dt r i (t) / Ylii&dt)=R ^- ^ n ^ ne °ther hand, it is obvious that any 
solution of Eqs. (2) is also be that of Eqs. (P) . In a mathematical language, Eqs. (P) and Eqs(2) 
have their respective solution sets: {solutions of Eqs. (2)} and {solutions of Eqs.(P)}, and 
there strictly exists a relation between these two sets: {solutions of Eqs. (2)} = {solutions of 
Eqs.(P)}. Namely, starting from the starting model equations of particle simulation scheme, 
we could find that there exists a closed equation set of Ufi, E and B, i.e., Eqs. (2. 0-4). 

b) for Vlasov-Maxwell theory. 

We could rewrite Vlasov equation (VE), for example Eq.(V.5), as 
= d t f + v-Vf-[E + vxB]-d p f. 

= [9t{f ~ fmono) + V ■ V (/ - f mono ) — [E + V X B] ■ d p (f — fmono)} 
+ (V-Ufi)- V fmono - (V-Ufi)- d p f mono 

+ [d t fmono +Ufi- V fmono ~ [E + Ufi X B] ■ d p f m ono] ■ (fV.l) 

Any distribution function / has two independent characteristic parameters: the variance 
and the mean. Here, the mean of / is represented by Ufi = J vfd 3 pj J fd 3 p). For any dis- 
tribution /, we could express it as / = nS (v — Ufi) + a 6 (v — Ufi) + Yli>i a i ( v ~ u fiY 
(where n = J fd 3 p, Ufi — J vfd 3 p/ J fd 3 p, a t are independent of v, a depends on 
all coefficients a^i through two relations, J a 6 (v — Ufi) + X)i>i a i ( v ~ u fiT d 3 p = 
and J v uo6 (v — Ufi) + Yli>i a i i v ~ u fiY d 3 p = 0, i.e. ao is a function of all a^i, 



). Here, a given pair of (n,Ufi) could correspond to multiple possi- 
ble distribution modes over w-space. This fact determines these two relations). Substituting 
this expression into VE and comparing the coefficients of (v — Ufi) l -term, we could find that 
there exists following equation for f mono — nS (v — Ufi) + a 6 (v — Uf t ) (because of the fact 
that VE is valid at any v- value.) 

= 8 t fmono +Ufl- Vf mono — [E + Ufl X B] • d p f mono , (fV.2) 

which could directly lead to (here, as stressed latter, p (ufi) = T n 2 ) 

v 1 ~ u n 

= d t [p (u fl )} + u fl ■ V r [p (u fl )} -[E + u fl xB] (fV.3) 

according to standard procedure, i.e., two relations, J [right-side terms in Eq.(fV.2)] d 3 p = 
and J p * [right-side terms in Eq.(fV.2)] d 3 p = 0, will lead to Eq.(fV.3), which is equivalent 
to Eq.(2.0) 

Indeed, because Eq.(fV.l) is of a binary-function type general form: = 
function(varl,var2), where varl and var2 are independent variables, like deriving 
Eq.(fp.2.b) from Eq.(fp.l), we could derive Eq.(fV.2) similarly (where var2 = v — Ufi) 

= J [right-side terms of Eq.(fV.l)] * 6 (var2) dvar2 
= right-side terms in Eq.(fV.2). (fV-2) 

c) for fluid theory 

It is well-known that fluid theory is indeed a derivant of V-M theory. All equations in 
fluid theory, except the assumed thermodynamical state equation, could be derived from 
5 equations in V-M theory according to standard procedure. For example, two relations, 
J [right-side terms in Eq.(fV.l)] d 3 p = and J p * [right-side terms in Eq.(fV.l)] d 3 p = 0, 
will lead to Eq.(F.5). 

Because the velocity v is a nonlinear function of the momentum p (i.e., v = J- — ) 

V i+p 2 

and vice versa, we should note that the statistic average value J pfd 3 p/ J fd 3 p (i.e. fluid 



momentum pfi) is usually not equal to the momentum corresponded by the statistic average 
value J vfd 3 p/ J fd 3 p (or fluid velocity Ufi), i.e., ^ p (u^) (where p(ufi) refers to the 
value of function p (variable) at variable = Ufi), if the distribution / is not a Dirac function 
of p (i.e., / has a thermal spread over p-space). Only at zero temperature case, there is 
Pfi = p(v,fi). (Strictly speaking, if / is a symmetric function of p, there will be pfi = 
p(ufi) — 0, Uj\ = and thermal pressure^ 0. But this special case corresponds to E = 
and E + Ufi x B = 0. A non-zero thermal pressure will drive p^ differing from according 
to Eq.(F.5). Once pfi ^ 0, there will be pfi ^ p (w/z) because / has an asymmetric thermal 
spread over p-space) . 

After noting the difference between pfi and p(ufi), a scrupulous reader will also note that 
it is Ufi, rather than pfi, that appears in Maxwell equations (Meqs). More important, he 
might consider whether or not there is necessity to introduce an assumed thermodynamical 
state equation. This is because two relations, J [right-side terms in Eq.(fV.2)] d 3 p = and 
J p * [right-side terms in Eq.(fV.2)] d 3 p = 0, will lead to Eq.(fV.3). Substracting Eq.(F.5) 
and Eq.(fV.3), an equation about thermal pressure/density will naturally appear. 

Therefore, we could strictly derive fluid theory from V-M theory according to a standard 
procedure without introducing any assumption. This makes fluid theory becoming really 
basic tool whose reliability is equal to those of other tools. In other words, the assumed 
thermodynamical state equation, Eq.(F.6), is replaced by Eq.(F.5) — Eq.(fV.3). Moreover, 
even if starting from particle simulation, we could still find that for all subindex i meeting 
i 6 rj (t) = R, summing corresponding Eq.(fp.2.a) will also lead to Eq.(F.5) — Eq.(fV.3). 
This also suggests that particle simulation and Vlasov-Maxwell theory completely agree with 
each other. 

By now, we have displayed in details how to obtain a closed equation set oiuji, E and B 
from three basic tools. In short, no matter which one of three basic tools is chosen by people 
when investigating plasma physics, E and B, obey a fixed fluid equation set, Eqs. (2.0-4). 
Indeed, these different basic methods are equivalent if they are in their respective strict 
forms. There is no reason to think that any method is better than others. 
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